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Abstract 

Three  consecutive  modes  of  flame  propagation  were  observed  over  a  bed  of  nano-aluminum  burning 
with  a  counter-flowing  oxidizer  of  20%  oxygen  and  80%  argon  by  volume,  each  displaying  significantly  dif¬ 
ferent  characteristics.  The  first  mode  of  propagation  was  the  focus  of  this  study  and  was  examined  within 
the  critical  Rayleigh  and  Peclet  number  regime  where  three-dimensional  buoyancy  effects  were  hindered 
and  the  fingering  thermal-diffusive  instability  occurred.  Fingering  flame  spread  was  observed  and  charac¬ 
terized  for  various  Peclet  numbers,  top  plate  heights  and  particle  sizes  to  gain  a  better  understanding  of  the 
reaction  mechanism  associated  with  the  combustion  of  nano-particles  in  close  contact.  Results  indicate  that 
the  first  mode  of  flame  propagation  over  a  bed  of  nano-aluminum  has  spread  rates  an  order  of  magnitude 
greater  than  that  of  cellulose  fuels.  However,  similar  trends  occur  when  varying  the  Peclet  number  and  the 
height  of  the  top  plate.  Furthermore,  faster  propagation  speeds  occur  with  smaller  particles  because  of 
their  increased  specific  surface  area.  The  widths  of  the  fingers  grow  and  more  of  the  surface  is  burned  with 
increasing  particle  size  due  to  the  longer  time  scale  available  for  lateral  growth. 

©  2006  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 

Keywords:  Flame  spread;  Nano-aluminum  powder;  Thermal-diffusive  instability;  Peclet  number;  Fingering  combustion 


1.  Introduction 

The  high  heat  of  reaction  of  aluminum  formu¬ 
lations  makes  it  a  performance  enhancer  in  solid 
rocket  propellants.  With  nano-sized  aluminum 
(nAl)  particles,  the  specific  surface  area  increases 
creating  easier  ignition  and  increased  burn  rates 
[1].  Nano-particles  can  also  be  mixed  with  metal 
oxide  powders  to  create  highly  energetic  metasta¬ 
ble  intermolecular  composites  (MICs)  for  a  varie- 


^  Supplementary  data  for  this  article  can  be  accessed 
online.  See  Appendix  A. 
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ty  of  applications.  MICs  can  have  burn  rates  of  up 
to  1  km/s  [2,3]. 

From  a  single  particle  [4-6],  to  an  ingredient  in 
a  composite  propellant  [7,8],  aluminum  has  been 
both  modeled  and  studied  experimentally.  How¬ 
ever,  in  order  to  gain  a  better  understanding  of 
its  combustion  properties,  aluminum  particles 
reacting  in  close  contact  is  an  area  that  needs  fur¬ 
ther  investigation.  When  a  pile  of  nAl  powder  is 
ignited  in  the  open  atmosphere,  a  flame  front 
spreads  across  the  surface.  In  contrast,  ignition 
and  flame  spread  cannot  easily  be  obtained  with 
micron-sized  aluminum  powder.  This  brings  up 
fundamental  questions  as  to  how  a  reaction  wave 
spreads  across  nAl  powder.  Past  flame  spread 
studies  of  solid  fuels  have  provided  a  controlled 
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manner  in  which  to  study  the  roles  of  fuel  and  oxi¬ 
dizer  properties  during  a  combustion  event. 

Flame  spread  over  thermally  thin  and  thick 
fuels  with  a  counter-flowing  oxidizer  have  been 
modeled  [9-12]  and  studied  experimentally  [13- 
15].  Experiments  by  Fernandez-Pello  et  al  [13] 
show  that,  at  low  oxidizer  flow  velocities,  the 
flame-spread  rate  over  thin  and  thick  fuels 
becomes  practically  independent  of  the  gas  veloc¬ 
ity.  It  is  hypothesized  that  at  these  low  velocities, 
buoyancy-driven,  three-dimensional  hydrody¬ 
namics  of  the  oxidizer  start  playing  a  role  in  feed¬ 
ing  the  flame. 

Investigations  of  flame  spread  in  microgravity 
with  low  oxidizer  velocities  have  shown  that  the 
spread  propagation  velocity  will  continue  to 
decrease  while  the  continuous  front  begins  to  sep¬ 
arate  into  “fingers”.  These  fingers  are  a  result  of  a 
thermal-diffusive  instability  present  at  these  low 
flows  [16,17].  Zik  and  Moses  [18]  were  able  to 
reproduce  these  instabilities  in  gravity  by  intro¬ 
ducing  a  top  boundary  to  the  flow.  This  top  plate, 
when  brought  to  a  critical  height,  or  critical  Ray¬ 
leigh  number  ( Ra ),  was  able  to  prevent  buoyancy 
flows  and  thus,  create  fingering  instabilities  in  a 
certain  Peclet  number  ( Pe )  regime.  Below  a  critical 
Ra  and  Pe ,  the  dimensions  and  spacing  of  the  fin¬ 
gers  could  be  well  controlled.  It  was  concluded 
that  the  distance  between  the  fingers  was  deter¬ 
mined  by  the  availability  of  oxygen,  while  the 
width  of  the  fingers  was  determined  by  the  heat 
losses  in  the  system.  The  propagation  speed  of 
the  fingers  depended  on  the  velocity  of  the  oxidiz¬ 
er  and  the  dimensions  of  the  fingers. 

The  objective  for  this  paper  is  to  examine  the 
fingering  flame  spread  of  a  bed  of  nAl  powder. 
Utilizing  nAl  powder  allows  for  investigations  of 
other  significant  variables,  such  as  particle  size 
(dp)  and  packing  density  (ppack)  that  cannot  be 
obtained  with  a  single  piece  of  solid  fuel  such  as 
filter  paper.  Moreover,  the  use  of  fingering  flame 
spread  allows  for  a  well-controlled  experiment  to 
give  insight  into  the  combustion  of  nano-alumi¬ 
num  particles  in  contact  with  each  other. 


2.  Experiment 

The  experiment  consists  of  a  copper  bottom 
plate,  copper  gas  diffuser  (5  in  Fig.  1)  and  an 
aluminum  bracket  which  holds  the  top  plate. 
For  visualization,  a  quartz  window  is  used  as 
the  top  plate  (1  in  Fig.  1).  To  prevent  lateral  flow 
over  the  bed,  copper  shims  were  placed  along  the 
sides  (2  in  Fig.  1).  The  top  plate  rests  on  the  shims, 
thus,  the  height  of  the  shims  also  determines  the 
height  of  the  top  plate  ( h ).  The  aluminum  sits  in 
the  bed  of  the  bottom  plate  (3  in  Fig.  1)  that  has 
dimensions  3.81  x  6.35  x  0.32  cm.  Bed  shims  with 
thickness  of  0.08  cm  were  also  created  to  adjust 
the  thickness  of  the  bed  (fbed).  The  gas  diffuser 


Fig.  1.  Schematic  of  the  nAl  bed  setup.  1.  quartz  top 
plate,  2.  copper  side  shims,  3.  bed  for  nAl,  4.  holes  for 
ignition  wire,  5.  copper  gas  diffuser. 


has  oxidizer  inlets  on  both  ends  and  a  row  of 
twenty-four  0.08  cm  diameter  holes  on  the  exit 
face  used  to  distribute  the  flow  across  the  width 
of  the  channel  and  assure  uniform  flow  develop¬ 
ment.  A  nichrome  wire  is  placed  through  holes 
at  the  opposite  end  of  the  bottom  plate  (4  in 
Fig.  1)  and  resistively  heated  by  a  DC  power  sup¬ 
ply  to  serve  as  an  ignition  source. 

To  capture  the  experiment,  a  Pulnix  TMC- 
6700-CF  camera  was  used  at  30  frames  per  second 
with  a  shutter  time  of  1/30  of  a  second.  Each 
frame  was  scanned  to  find  the  v-location  of  pixels 
with  intensities  above  10%  of  the  maximum. 
These  values  were  recorded  and  averaged  to 
acquire  an  average  v-location  for  each  frame. 
The  average  v-location  and  the  time  between  each 
frame  were  used  to  calculate  an  average  velocity 
of  the  front  (rf).  Depending  on  the  experimental 
condition,  an  appropriate  aperture  was  chosen 
to  capture  enough  light  without  saturation. 

The  oxidizer  was  a  gas  mixture  of  20%  oxygen 
and  80%  argon  by  volume.  To  control  the  flow, 
two  mass  flow-controllers  (Hastings  HFC-202) 
were  used  (one  for  the  oxygen  and  one  for  the 
argon).  The  oxygen  fraction  and  flow  velocity 
were  controlled  by  a  computer  data  acquisition 
system  (Measurement  Computing  USB-1208FS 
DAQ  Board  and  National  Instruments  Fab  view 
7.1).  Fluid  dynamics  calculations  (Fluent  6.2) 
were  performed  on  the  extreme  cases  to  find  an 
adequate  distance  for  the  flow  to  become  fully 
developed.  The  distance  from  the  exit  of  the 
gas  diffuser  to  the  beginning  of  the  nAl  particle 
bed  was  found  to  be  2.5  cm  to  assure  fully  devel¬ 
oped  flow.  The  average  oxidizer  velocity  (uox) 
was  determined  by  the  channel  cross-sectional 
area  and  the  oxidizer  flow  rate.  For  some  experi¬ 
ments,  the  entire  exhaust  gas  exiting  from  the  test 
section  was  sampled  and  analyzed  by  gas  chroma¬ 
tography  (Agilent  Micro  GC3000)  for  oxygen 
concentration. 

A  pre-measured  mass  of  nAl  powder  was 
placed  in  the  bed  and  distributed  uniformly.  To 
create  a  smooth  top  surface,  a  glass  surface  was 
used  to  lightly  press  the  powder  into  the  bed. 
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Four  different  sized  particles  were  used  having 
diameters  of  38  (Technology  Materials  Develop¬ 
ment),  50,  80,  and  120  nm  (latter  three  from 
Nanotechnologies,  Inc.).  The  oxide  passivation 
layers  are  3.1,  2.1,  1.9,  and  1.8  nm  and  the  pyc¬ 
nometer  densities  are  3.2,  3.0,  3.1,  and  2.9  g/cm3, 
respectively.  An  SEM  micrograph  of  80  nm  nAl 
in  Fig.  2  shows  the  approximate  spherical  shape 
and  the  uniformity  of  the  particles. 

Experiments  were  conducted  specifically  to 
study  the  effects  of  Pe ,  h,  and  dp  on  Vf,  the  average 
width  of  fingers  (wf),  and  the  average  distance 
between  fingers  (wb).  The  Peclet  number  was 
defined  as  v0Xh/D02 ,  where  D0l  is  the  molecular 
diffusion  coefficient  of  oxygen  and  assumed  con¬ 
stant  throughout  at  a  value  of  0.25  cm2/s. 


3.  Results  and  discussion 

We  observe  three  modes  of  flame  propagation 
(Fig.  3)  that  occur  sequentially  throughout  each 
experiment:  1 .  counter-flow  surface  reaction  wave, 
2.  co-flow  bulk  reaction  wave,  and  3.  cellular 
flames  randomly  propagating  through  the  materi¬ 
al  somewhat  similar  to  what  was  observed  by 
Zhang  et  al  [14]  with  filter  paper  in  low  Lewis 
number  atmospheres  near  the  extinction  limit. 
Son  et  al  [19]  observed  the  bottom  of  a  pile 
of  nAl  during  the  first  wave  to  confirm  this  was 
only  a  surface  burn.  This  study  focused  on  the 
first  wave. 

To  demonstrate  the  effect  of  Pe  and  A,  38  nm 
particles  were  packed  at  equal  ppack  of  0-22  g/ 
cm3  in  beds  with  a  thickness  (/bed)  of  0.32  cm. 
The  effect  of  dp  was  investigated  at  two  different 
Pe  values  and  h  =  2  mm  with  particle  diameters 
of  50,  80,  and  120  nm.  A  ppack  of  0.34  g/cm3  was 
used  for  these  experiments.  Due  to  the  inherent 
pouring  density  of  the  Nanotechnologies  Inc.  par¬ 
ticles,  a  larger  mass  was  needed  to  completely  fill 
the  bed.  Because  of  a  limited  amount  of  these  par¬ 
ticles  available,  thed  was  reduced  to  0.08  cm  for 
these  experiments.  Experiments  showed  this  did 
not  affect  the  first  wave. 

Figure  4a  shows  a  sequence  of  video  frames 
(A/  =  1/15  s).  The  flame’s  trail  loses  luminosity  a 


Fig.  2.  SEM  micrographs  of  80  nm  aluminum  particles 
manufactured  by  Nanotechnologies,  Inc. 


Fig.  3.  Three  modes  of  combustion  for  a  bed  of  38  nm 
nAl. 


Fig.  4.  (a)  Sequence  of  frames  each  1/15  s  apart  (b)  all 
frames  integrated  through  the  run  to  show  the  history  of 
the  flame  spread  for  a  38  nm  nAl  bed. 


Fig.  5.  Integrated  flame  histories  of  beds  of  38  nm  nAl. 
Refer  to  Fig.  4b  for  a  Peclet  number  (Pe)  of  3  and  top 
plate  height  ( h )  of  2  mm  case. 


certain  distance  behind  the  front.  This  trail  length 
is  the  order  of  1  mm  long.  Integration  of  all  of  the 
frames  along  the  whole  length  of  the  bed  (Fig.  4b) 
shows  the  history  of  the  flame  spread.  Along  with 
Fig.  4b,  Fig.  5  shows  flame  histories  for  a  bed  of 
38  nm  particles  for  Pe  =  3,  9,  and  16  at  h  =  2 
and  for  h  =  2,  6,  and  10  mm  at  Pe  =  3. 

3.1.  Conservation  of  mass 

Analysis  of  the  exhaust  gases  by  gas  chroma¬ 
tography  at  the  exit  during  the  first  wave  showed 
an  oxygen  concentration  of  0.8%,  a  96%  decrease 
from  the  original  amount.  It  is  reasonable  to 
assume  in  the  analysis  that  all  of  the  oxygen  is 
consumed  by  the  front  and  a  global  mass  conser¬ 
vation  can  be  written  as 
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—  /4toicA)2,  (1) 

where  r Ai  is  the  mass  consumption  rate  of  alumi¬ 
num  per  unit  area,  r0l  is  the  mass  consumption 
rate  of  oxygen  per  unit  area  and  psto ic  is  the  stoi¬ 
chiometric  coefficient. 

When  varying  the  height  of  the  top  plate,  Pe  is 
kept  constant.  Therefore,  the  mass  flow  rate  per 
unit  area  of  oxidizer  (mJJx)  and,  thus,  r0l  are  con¬ 
stant.  Since  r o2  is  constant,  r Ai  must  be  equal  for 
each  case  according  to  Eq.  (1). 

For  the  Pe  =  3,  h  =  6  mm  and  Pe  =  3,  h  = 
10  mm  cases  the  first  40%  of  the  burn  is  neglected. 
Notice  from  Fig.  5  that  this  area  contains  notice¬ 
ably  more  burned  surface  area  than  the  latter 
60%.  At  these  Peclet  numbers,  the  vox  is  relatively 
slow.  The  slow  velocity  could  possibly  allow  a 
small  amount  of  air  to  diffuse  through  the  open 
end  at  the  exit  of  the  apparatus.  An  order  of  mag¬ 
nitude  analysis  shows  that  this  could  only  have  a 
small  effect,  but  preliminary  experiments  with  no 
oxidizer  flow  also  show  this  phenomenon.  More¬ 
over,  since  the  oxygen  is  being  consumed,  the 
gas  pressure  will  decrease  after  the  front  creating 
a  pressure  gradient  that  could  induce  flow  in  from 
the  ambient  at  a  higher  rate  than  estimated  by  dif¬ 
fusion.  These  observations  could  also  be  an  initial 
transient  for  the  instability  that  is  exaggerated  by 
the  high  top  plate  heights  and  low  advective  flows. 

The  r Ai  is  calculated  using 

TAi  =/bPAibed^f%,  (2) 

where  pAibed  is  the  packing  density  of  just  the 
aluminum  in  the  bed  (ppack  without  oxide  layer), 
/b  is  the  fraction  of  surface  area  burned,  and  <5b 
is  the  non-dimensional  depth  of  burned  aluminum 
scaled  by  a  particle  diameter  (<5b  =  dp/dref  where 
dre f  is  a  reference  depth  of  50  nm,  a  typical  nAl 
particle  size).  A  non-dimensional  depth  is  used 
because  the  depth  of  the  burn  is  not  known  for 
this  experiment.  It  is  assumed,  however,  that  the 
particles  lined  up  directly  on  top  of  each  other 
and  the  depth  of  the  burn  scaled  with  the  diameter 
of  the  particle.  The  TA i  was  found  to  be  0.42,  0.43, 
and  0.49  g/cm2-s  for  h  of  2,  6,  and  10  mm,  respec¬ 
tively.  The  TA1  for  h  =  10  is  slightly  higher  due  to 
the  slowest  oxidizer  velocity  further  exaggerating 
the  initial  transient  discussed  above. 

When  varying  Pe ,  h  is  kept  constant.  As  Pe 
increases,  TA1  should  increase  accordingly.  Equa¬ 
tion  ( 1 )  can  be  simplified  to  /brf  =  Zvox ,  where  Z 
is  a  constant  equal  to  /WPo^ox/Paw  aox  is 
the  non-dimensional  height  of  the  oxidizer  flow 
channel  (Aox  =  /z//zref  where  hre f  is  a  reference 
height  of  2  mm,  the  lowest  h  used  in  the  experi¬ 
ment),  and  po,  is  the  density  of  oxygen.  Therefore, 
Vf  multiplied  by  /b  should  vary  linearly  with  the 
velocity  of  the  oxidizer  as  shown  in  Fig.  6. 

The  values  of  pA1,  are  found  to  be  0.29,  0.26, 
and  0.24  g/cm3  for  dp  of  50,  80,  and  120  nm, 
respectively.  The  TA1  can  be  found  for  each  case 


Fig.  6.  Relation  between  the  flame  propagation  velocity 
(ttf)  and  the  oxidizer  velocity  (vox)  for  varying  the  Peclet 
number  (Pe). 


at  the  same  Pe  and  compared  using  Eq.  2.  For 
Pe  =  3,  Jai  =  0.06,  0.07,  and  0.06  g/cm2-s  for 
dp  =  50,  80,  and  120  nm,  respectively.  For 
Pe  =  5,  Jai  =  0.16,  0.11,  and  0.13  g/cm2-s  for 
dp  =  50,  80,  and  120  nm,  respectively. 

3.2.  Effect  of  top  plate  height 

Figure  7  shows  the  trends  for  vvf,  wb,  and /b  at  a 
constant  value  of  Pe  =  3  with  a  varying  h  value. 
Since  the  Pe ,  and  thus,  the  mass  flow  rate  per  unit 
area  of  oxygen  (wq  )  are  constant,  the  /b  remains 
nearly  constant  as  well.  There  is  a  significant 
increase  in  vvf  with  an  increase  in  h,  while  vvb 
remains  nearly  constant.  From  the  conservation 
of  mass  analysis  above,  v f  must  stay  constant  as 
is  shown  in  Fig.  8.  Overall,  when  h  increases,  vox 
must  decrease  to  keep  Pe  constant.  With  a  slower 
vox  and  a  higher  h  the  flame  front  losses  some  abil¬ 
ity  to  conduct  heat  out  through  the  top  plate  due 
to  a  decreasing  convective  flow  and  an  increase  in 
the  length  scale  for  conduction.  This  loss  of  ability 
for  the  front  to  conduct  heat  out  through  the  top 
results  in  an  increase  in  wf .  To  compensate  for  the 


h  (mm) 

Fig.  7.  Trends  for  the  width  of  the  fingers  (wf),  width 
between  the  fingers  (wb),  and  the  fraction  burned  (fh)  at 
Peclet  number  (Pe)  of  3  for  38  nm  nAl  bed. 
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Pe 


Fig.  8.  Trend  for  flame  propagation  velocity  (uf)  with 
varying  top  plate  height  ( h )  or  Peclet  number  ( Pe )  for 
38  nm  nAl  bed;  insert  shows  the  steady  nature  of  the 
front. 

increase  in  wf,  fewer  fingers  are  present  with  prac¬ 
tically  constant  values  of  wb,  /b,  and  Vf.  The  aver¬ 
age  flame  front  velocity  is  constant  throughout  the 
run  as  shown  in  the  insert  in  Fig.  8.  Each  finger 
will  remain  at  a  velocity  nearly  equal  to  the  overall 
average  Vf  if  it  is  undisturbed.  However,  if  a  finger 
propagates  in  front  of  another  it  will  “screen”  the 
finger  from  the  oxygen  and  extinguish  it. 

3.3.  Effect  of  Peclet  number 

Figure  9  shows  that  /b  increases  as  Pe  increas¬ 
es,  which  results  from  the  increase  in  /Wq  of  the 
oxidizer  stream.  For  the  Pe  =  3  and  Pe  =  9  cases, 
Wf  was  found  nearly  constant  increasing  only 
slightly  from  3  and  5  mm,  respectively.  The  work 
of  Zik  et  al.  [18]  showed  this  trend  as  well  with 
filter  paper  for  a  constant  h.  Since  these  conditions 
have  nearly  constant  values  of  Wf,  wb  must 
decrease  when  /b  increases  as  is  shown  in  Fig.  9. 
At  Pe  =  16,  a  nearly  flat  front  is  observed  forcing 


Fig.  9.  Trends  for  the  fraction  burned  (fb)  and  the  width 
between  the  fingers  (wb)  at  a  top  plate  height  ( h )  of  2  mm 
for  38  nm  nAl  bed. 


Fig.  10.  Integrated  flame  spread  histories  for  different 
particle  diameters  (<ip)  at  Peclet  numbers  (Pe)  of  3  and  5 
and  a  top  plate  height  ( h )  of  2  mm. 


Wf  to  be  the  width  of  the  bed,  wb  to  be  zero,  and 
the  total  surface  area  to  be  consumed.  For  this 
condition,  the  critical  Pe ,  which  is  a  threshold  val¬ 
ue  where  the  fingering  instability  ceases  to  occur, 
has  been  exceeded.  Figure  8  shows  that  Vf  increas¬ 
es  significantly  with  increasing  Pe.  This  increase  is 
directly  related  to  the  conservation  of  mass  analy¬ 
sis  for  increasing  Pe. 

3.4.  Effect  of  particle  size 

Figure  10  shows  integrated  histories  of  the 
flame  spread  for  50,  80,  and  120  nm  at  Pe  =  3 
and  Pe  =  5.  Since  smaller  aluminum  has  a  higher 
specific  surface  area,  the  burn  time  for  particles 
will  decrease  with  decreasing  diameter,  which 
should  lead  to  faster  front  velocities  for  beds  of 


—  A  -  -  38  nm  H 80  nm 

- © 50  nm  0  120  nm 


Pe 


Fig.  11.  Trend  for  flame  propagation  velocity  (rf)  with 
varying  Peclet  number  (Pe)  for  different  particle  diam¬ 
eters  (dv)  at  a  top  plate  height  (h)  of  2  mm. 
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smaller  particles.  These  velocities  are  shown  in 
Fig.ll.  The  data  from  the  38  nm  cases  further 
demonstrate  this  trend  even  though  these  particles 
have  a  larger  oxide  passivation  layer,  lower  ppack, 
and  a  larger  /bed. 

Changing  the  particle  size  will  also  have  an 
effect  on  Wf  and  /b.  Since  a  bed  of  smaller  particles 
will  propagate  faster  due  to  the  higher  specific  sur¬ 
face  area,  the  /b  must  decrease  according  to  Eq. 
(2).  This  trend  is  shown  in  Fig.  12.  For  the  larger 
particles,  the  slower  front  velocities  allow  more 
time  for  the  lateral  diffusion  of  the  flame  causing 
wider  fingers.  Furthermore,  because  there  are  few¬ 
er  particles  per  bed  for  the  same  ppack,  larger  par¬ 
ticle  sizes  will  have  less  contact  resistance.  This 
will  lead  to  a  higher  thermal  conductivity,  which 
could  aid  in  the  widening  of  the  fingers. 

3.5.  Other  observations  and  future  work 

Notice  in  Fig.  11,  Vf  becomes  independent  of 
Pe  (or  rox)  at  Pe  >  9  for  the  slower  burning  parti¬ 
cles.  At  these  higher  values  of  Pe  the  whole  sur¬ 
face  of  the  nAl  bed  is  burned  (fh  — »  1,  and  wh 
— >  0).  Thus,  the  mass  of  oxygen  available  begins 
to  exceed  the  stoichiometric  value  and  not  all  of 
it  is  consumed  by  the  front.  Therefore,  the  front 
is  “saturated”  and  has  reached  a  maximum  v{. 
This  independence  of  oxidizer  velocity  indicates 
that  the  nAl  bed  is  thermally  thin  according  to 
deRis’  [9]  model  of  opposed  flow  flame  spread. 
However,  Vf  is  also  independent  of  /bed?  as 
described  above,  indicating  that  the  bed  is  ther¬ 
mally  thick  [9].  The  explanation  for  this  contradic¬ 
tion  is  that  the  whole  bed  is  not  participating  in 
the  first  wave.  It  can  be  considered  a  thermally 
thin  flame  spread  resting  on  top  of  a  non-partici¬ 
pating  bed  of  fuel. 


■IS -  fb  at  Pe=5  - 0 -  wf  at  Pe=5 

0  fb  at  Pe=3  - © -  wfat  Pe=3 


Fig.  12.  Trends  for  the  width  of  the  fingers  (Wf)  and 
fraction  burned  iff)  with  varying  particle  diameters  (dp) 
or  surface  to  volume  ratios  (S/V)  at  two  different  Peclet 
numbers  for  a  top  plate  height  (h)  of  2  mm. 


The  integrated  histories  of  each  particle  size 
displayed  distinct  fingering  patterns  at  certain 
Peclet  numbers.  Figure  10  shows  the  experiments 
with  Pe  =  3  having  similar  structures.  However, 
with  Pe  =  5,  the  fingers  look  significantly  different 
depending  on  dv.  The  80  nm  particle  bed  has  fin¬ 
gers  that  have  a  feathery  look,  for  example.  The 
branches  are  wider  than  the  50  nm  particle  bed 
and  the  edges  have  a  much  smoother  intensity  gra¬ 
dient.  The  120  nm  particle  bed  has  fingers  that  are 
not  well  defined  and  the  edges  are  extremely 
jagged. 

During  some  runs,  the  front  seems  to  jump  for¬ 
ward  as  shown  in  Fig.  13.  This  could  be  an  exam¬ 
ple  of  an  agglomeration  of  particles  exploding  and 
propelling  the  hot  products  forward.  Particle 
explosion  was  observed  by  Yetter  et  al.  [4]  for  sin¬ 
gle  particles.  More  frequent  explosions  of  this 
magnitude  were  noticed  with  smaller  particles  at 
higher  propagation  velocities,  although  this  mode 
of  propagation  is  the  exception. 

An  experiment  was  performed  with  no  top 
plate  and  slow  vox.  With  no  top  plate,  three-di¬ 
mensional  buoyancy  effects  are  not  suppressed 
and  the  existence  of  fingers  would  not  be  antici¬ 
pated  based  on  previous  experiments  with  filter 
paper  combustion.  However,  fingers  are  present 
even  without  a  top  plate,  which  could  be  due  to 
a  faster  Vf  for  nAl  than  for  filter  paper  allowing 
less  time  for  buoyancy  flow  to  develop.  There 
was,  however,  more  surface  area  burned,  more 
fingers  and  a  faster  Vf  than  with  the  top  plate  indi¬ 
cating  that  the  restricted  buoyant  flow  from  the 
top  plate  does  play  a  significant  role  in  the  finger¬ 
ing  behavior.  The  induced  buoyant  flow  causes  Vf 
and  the  finger  dimensions  to  be  similar  to  a  run 
with  Pe  =  9. 

After  the  first  and  second  waves,  no  white  alu¬ 
mina  powder  was  observed,  which  is  what  would 
be  expected  if  stoichiometric  A1203  is  formed. 
Only  after  the  third  mode  were  trails  of  white 
noticed.  It  is  hypothesized  that  the  first  wave  only 
thickens  the  already  present  A1203  layer  a  small 
amount  (possibly  changing  its  polymorph  struc¬ 
ture)  to  remain  optically  thin  and,  therefore,  the 
particles  do  not  appear  white  in  color,  but  are 
slightly  discolored  from  the  unburned  bed.  The 
second  wave  could  leave  the  particles  optically 
thin  as  well,  or  since  it  is  diffusion  limited,  there 
could  be  non-stoichiometric  products  formed. 
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Fig.  13.  Example  of  particle(s)  explosion  in  80  nm  nAl 
bed. 
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4.  Conclusions 

Flame  spread  was  examined  over  a  bed  of  nAl 
particles  with  20%  02  and  80%  Ar,  by  volume,  as 
the  oxidizer.  Similar  to  the  experiment  of  Zik  et  al. 
[18],  a  top  plate  was  used  to  inhibit  three-dimen¬ 
sional  buoyancy  effects  and  create  a  fingering 
instability  that  was  only  observed  in  flame  spread 
experiments  in  microgravity.  This  gave  insight 
into  the  nature  of  nano-aluminum  particles  burn¬ 
ing  in  close  contact  with  each  other. 

Specifically,  we  found: 

1.  Three  flame  propagation  modes  are  observed 
over  a  bed  of  nAl  with  an  opposed  oxidizer 
flow,  including  a  discrete  cellular  flame  mode. 

2.  The  first  mode  of  flame  propagation  across  a 
bed  of  nano-aluminum  particles  was  investi¬ 
gated  in  this  study  and  has  similar  characteris¬ 
tics  to  a  continuous  solid  fuel,  although  the 
spread  rate  is  much  more  rapid. 

3.  For  an  increasing  Peclet  number,  the  width 
between  the  fingers  decreases,  the  percent 
burned  of  surface  area  increases,  and  the  prop¬ 
agation  velocity  increases. 

4.  For  an  increasing  top  plate  height,  the  width 
between  the  fingers  remains  constant,  the 
width  of  the  fingers  increases,  and  the  propaga¬ 
tion  velocity  remains  constant. 

5.  Increasing  particle  size  decreases  the  propaga¬ 
tion  velocity  and  widens  the  fingers.  The 
decreased  propagation  velocity  is  due  to  the 
lower  specific  surface  area  of  the  particles. 
Because  of  this  decreased  propagation  velocity, 
the  fingers  have  more  time  to  diffuse  laterally 
creating  wider  fingers.  Furthermore,  the  higher 
thermal  conductivity  of  larger  particle  beds 
likely  aids  in  the  lateral  diffusion  of  the  flame. 

6.  The  flame  history,  particularly  at  higher  Peclet 
numbers,  showed  distinct  finger  shapes 
depending  on  the  size  of  the  particle. 

7.  Particles,  or  perhaps  groups  of  particles,  were 
observed  to  explode  causing  the  front  to  jump 
ahead  in  certain  areas. 

8.  Removing  the  top  plate  did  not  eliminate  the 
fingering  in  contrast  to  filter  paper  flame 
spread.  The  flame  spread  rate  and  finger 
dimensions  were  similar  to  that  of  a  run  at  a 
Peclet  number  of  9. 
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Comment 

Andrew  Sullivan,  Australian  National  University,  Aus¬ 
tralia.  Did  you  try  initial  flame  spread  in  the  same  direc¬ 
tion  as  oxidizer  (forward  flow)  and  would  you  expect  to 
see  different  behavior  regarding  fingering? 

Note:  fingering  resembles  dendritic  behavior  of  some 
cellular  automata  of  diffusion  limited  aggregation 
modeling. 

Reply.  I  did  not  try  concurrent  flame  spread,  how¬ 
ever,  Olson  et  al.  [1]  performed  this  in  a  similar 
experiment  with  paper  and  found  that  the  fingers 
moved  back  and  forth  laterally,  had  a  less  distinct 
shape,  burned  more  of  the  fuel  surface,  and  had  a 
lower  steady  state  velocity.  I  would  expect  this  to 


be  similar  with  nano-aluminum  because  it  is  attribut¬ 
ed  solely  to  hydrodynamic  effects  which  are  similar 
for  the  two  systems. 

Thank  you  for  the  comment  on  similarity  to  cellu¬ 
lar  automata  of  diffusion  limited  aggregation  model¬ 
ing.  In  [2],  they  relate  this  phenomenon  to 
electrochemical  deposition,  another  form  of  diffusion 
limited  growth  of  fingers. 
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